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Abstract 

We calculate ^/s dependence of Cronin momentum for p + A and A + A collisions 
in saturation model. We show that this dependence is consistent with expectation from 
formula which was obtained using simple dimentional consideration. This can be used 
to test validity of saturation model (and distinguish among its variants) and measure x 
dependence of saturation momentum from experimental data. 

1 Introduction 

It was shown in 1 that saturation model can explain Cronin like behavior. As many 
other models also explain this behavior we need some subtle prediction to distinguish one 
model from another. One prediction of this type is to calculate position of maximum in 
Cronin ratio. But there are three variable in Cronin effect to measure: momentum where 
Cronin ratio have maximum (let's call it Cronin momentum qc), value of maximum Rc 
and momentum where Cronin ratio equal to unity So why maximum position? We 
know that value of Cronin ratio in yl + ^ collisions have normalization uncertainty. And 
therefore variables Rc, Qu are not a good ones to make predictions. On the other side 
Cronin momentum does not depend on normalization and therefore is the best candidate 
we have. Let's consider for now central rapidity p + A collisions only. In saturation model 
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there is only one semihard scale which govers momentum dependence of differential cross- 
section ^^p- (here y rapidity and q transverce momentum of produced particles). It is 
saturation momentum Qg. As we have only one semihard scale then Cronin momentum 
can only depend on this scale. Using dimentional consideration the only equation which 
relates Qs and qc can be: 

qc = PQs (1) 

where f3 is some dimentionless constant. But as we know saturation momentun Qs is not 
a constant. It depends on Bjorken variable x which in this process defined by relation 

and as the only known scale is qc then instead of we'll have 

qc = PQs ( ^) (3) 



where /3i is another dimentionless constant. 

It is easy to define Qs{x). From geometric scaling effect for small x we have 

Ql{x)=Ay^Ql,(^^)' (4) 

where A = 0.3 is geometric scaling constant and QsO,xq some parameters those exact 
value we define from fact that for reaction with Au nucles at ^/s = 200 Get" we have 
= 1 — 2 Gev. It is easy to solve Q and write expression for qc 

qc = g^^3(2^Vs^ (5) 

or if we log both parts we'll have: 

\n{qc)=a + b\n{^s) (6) 

where a and b defined as: 

^ ln{A) + ln{ql) (7) 
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So to test saturation model prediction we should calculate Cronin momentum for different 
energies and check ®. There is however soft scale Aqcd and it is not obvious that it's 
existance does not change this formula. So we should check this dependence explicitly 



by numerical calculation. Let us consider Cronin ratio for p + A (from here we take Au 
nucleus) collision 



dyd-'q 

A- 



dyd''q 

As we stated before we suppose y = 0. In saturation model gluon production cross-section 
can be expressed as 

^^2^ = 7 ^ )'^^(^2,(g-/^) ), (9) 

where 4>a,p is unintegrated gluon distribution of nucleus and proton and xi,X2 defined by 
equation 

XI = -^e-y, X2 = -^ey (10) 

In leading logarifmic order we can rewrite @ in following form 
dfjP^ 2 a, 1 



{xiGp{xi,q^)(j)A{x2,q^) + X2GA{x2,q^)(t)p{xi,q^)^ , 



ill) 



d?q dy Cp q^ 

where xG{x,q'^) is gluon distributin function which can be expressed from (/)(x,fc^) by 
following relation 

xG(x,q'^) = r ct){x,k^)dk^ (12) 
Jk 

Using the same approximation we can write following approximate relation for cross-section 
of gluon production va. p + p collisions 

= ~c^~j^ {xiGp{xi,a^)(l)p{x2,q^) + X2Gp{x2,q^)(t)p{xi,q^)^ , (13) 
Let us suppose that in considered kinematical region unintegrated gluon distibution func- 



tion of proton (j)p{x,q^) does not depend on x and that (f)p{x,q^) = ^ ■ Then Cronin 
ratio can be expressed as: 



P A\ct>p{x2,q^) Gp{x2,q^) J 



or as we supposed y = then 

_ 1 [ (t>A{x,q'^) , GA{x,q'^) \ 

A[Mx,q^)^ Gp{x,q^)) ^''^ 

where x = -7= 

All we need now is the expression for unintegrated gluon distribution function. We 
consider three models for gluon distribution function: Kharzeev-Levin-Nardi proposed in 
P], MacLerran-Venugopalan proposed in 0J|S] and "dipole" model. 



2 Kharzeev-Levin-Nardi model 



In simplified form of this model unintegrated gluon distribution function (j){x,q'^) have 
following form: 



(16) 



where 4>\ normalization factor. 



Then for gluon distribution function G{x,q'^) we have: 

xGa{x, q^) = (IPa {q^ - ^qcd) , Q < Qs{x) 



(17) 



xGA{x,q^)=^''A\Qlix)ln 



Qlix) 



+ Q'l{x)-hlcD h9 > Qs{^ 



And for Cronin ratio we have (we use here approximate formula (|15jl ) 
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,q<Qs{x) (18) 



1 Q TT 

RpA = -^<t>A—pr- 
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,q>Qs{x) (19) 



If we look at (dni) we'll see that RpA here is non-decreasing function of momentum q so is 
not clean if there is any Cronin like behavior in this model. It should be mentioned however 
that x depends on q by means of relation x = (if we suppose y = 0) and as saturation 
momentum in low x region depend on x as Qlix) = A^/^QIq i^)^ then in reality ^ have 
maximun at some monentum qc (figure^ which value approximately defined by equation 



qc = Qs 



qc 



(20) 



and modified slightly by logarifmic terms in H19|l . Nevertheless we have formula similair to 



In(gc') = a + 61n(-v/i) 



(211 



where 6 = 0.1042 
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3 McLerran-Venugopalan model 

In McLerran-Venugopalan model expression for unintegrated gluon distribution function 
was finded in works [H] and can be written as 

ct>A{x,q') = / dhd\e~^^- 1 (1 - e-^'^'^'<^'^^^y% (22) 

Or if consider cilindrical nucleus 



</,^(x,g2) = f dVe-^- 1 (1 - e-^«^i-(iAA)/4), (23) 

or 



<^A{x,q^) = / drMq-r) 1 (1 - e^^^? '"(V^^V^), (24) 

as (27r)^ J r 

It is better however use the expression proposed in 8^ which relates unintegrated gluon 

distribution function in McLerran-Venugopalan model and the forward amplitude of sca- 

tering NG{r,x) of a gluon dipole of transverse size r and rapidity y = ln(l/x) on nucleus. 

When we can rewrite previous equation in the following form: 

Mx,q'') = I drUq.r) \ Ng{t,x), (25) 

where J^i^x) is Bessel function. It obvious that Nq depends on x. There are different ways 
to set this dependence. We can use Balitsky-Kovchegov equation to define x dependence 
of dipole scattering cross-section but as it is unsolved for now we choise more simple way. 
Let us define ad hoc that 

Nckr, X) = 1 - (.-'■''Qs{^?ln{XlrK)lA (26) 

i.e. all X dependence goes in definition of Qs(x). But we can not use (|26() directly as Ndr, x) 
have not very good behavior for large r( i.e. if r — > cxd then NG{r,x) becomes negative 
instead of unity). So we should regularize (|26() somehow. Let us regularize NQ{r,x) by 
following prescription: 

Ncir, x) = l- e^'Q4^')'('«('^A)-V(ln(rA))2+.2)+ln(roA))/8 (27) 

and set tq = and e < l(the final result does not depend on exact value of e, if it is not 
too large). It should be noted that result does not depend on regularization scheme and 
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we could regularize NQ{r,x) with something like this: 

Ncir, x) = l- e-^'-'Q^(^)''"(iAA)/4 ^ ^ ^ (28) 
NGir,x) = 1 - e--'Q4-)^in(iAoA)/4 > 

but this regularization is inconvenient in "dipole" model. 
Then for gluon distribution function G(x,g^) we have: 

xGA{x,q^) = ^j^2 j dr{qJ,{q-r)-AJ,{A-r))^NG{r,x)), (29) 

If we substitute this functions in Cronin ratio we'll have dependence which is presented 
in figure 121 Then we can calculate numericaly value of Cronin momentum qc for different 
energies. The result is presented in figure |1J The slope is b = 0.1323. It should be 
mentioned that even the line have different position they have almost the same slope as 
previous model and almost exactly equal one which was calculated in ((T)). 

4 "Dipole" model 

In "dipole" model we can relate unintegrated gluon distribution function with gluon dipole 
cross-section. It was done in work (HI |7] and expression for unintegrated gluon distribution 
function can be written as (we supposed that nucleus is cilindrical) 

Mx,q') = J d^e-'i-^' Vl NG{r,x)), (30) 

or 

M^. i) = / dr Mq ■ r) rV^. Ncir, x), (31) 

as {2TTy J 

For for gluon distribution function G{x,q'^) we'll have 

xGA{x,q^) = ^yf, 2 / dr{qJi{q-r)-AJi{A-r)) V^,NG{r,x), (32) 
as (27r)^ J 

As before it is easy to calculate numerically value of Cronin ratio. Result is presented on 
figure 01 Varying energy we calculate numerically Cronin momentum qc (result presented 
on figureEJ- And as before we have dependence In(gc') = a+61n(y^), with slope b = 0.1120 

5 A+A collisions 

Like in p + A collision in ^ + ^ collisions (we take only central rapidity region) there is 
only one semihard scale Qs- Ant therefore dependence of Cronin momentum qc must be 
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p + A 


A + A 


Model 






Kharzeev-Levin-Nardi 


0.1042 


0.1485 


McLerran-Venugopalan 


0.1323 


0.1383 


'Dipole' 


0.1120 


0.1244 



Table 1: Summary of slopes for different models in p + A and A + A collisions. 

govered by ©. We can apply all formulas above to this case, as we have for Cronin ratio 
following approximate relation similair to (|15)) 

Gp{x.p) (t)p{x,p) 

calculating numerically q dependence of Cronin ration for considered models (figures r5l6l7|) 
at differen energies we have same linear behavior for ln(qc) as before (figure |HJ and also 
have slopes consistent with (O. All data summarized in Table 1 (for 'dipole' model only 
points with ^/s > 500Gev was taked for slope calculation). 



6 Conclusion 

We calculate ^/s dependence of Cronin momentum in several models based on saturation 
and show that this dependence is consistent with simple formula based on geometric scaling 
only. This subtle prediction can test validity of saturation model. Even more. As slope 
values is slightly different we have posibility to distinguish among variants. But this requires 
more precice measurement of Cronin effect (at least in midle momentum region) that those 
we have today. Having this we can in turn measure dependence of saturation momentum 
Qsix) on X. 
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Figure 1: Cronin ratio for Kharzeev-Levin-Nardi gluon distribution function for p + A collisions 
y/s = 2006*61; (solid curve) and y/s = ITOOGev (dashed curve) 




Figure 2: Cronin ratio for McLerran-Venucopalan gluon distribution function for p + A collisions 
for = 200(761; (solid curve) and a/s = 1700G6V (dashed curve) 
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Figure 3: Cronin ratio for "dipole" gluon distribution function for p + A collisions for ^/s 
200Gev(soild curve) and ^/s = ITOOGef (dashed curve) 




Figure 4: Dependence of ln{qc){foT RpA)on ln{y/s) for different models: Kharzeev-Levin- 
Nardi(solid curve), McLerran-Venugopalan(dashed curve), "dipole" (dot-dashed curve) 
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Figure 5: Cronin ratio for Kharzeev-Levin-Nardi gluon distribution function for ^ + A collisions 
a/s = 2006*61; (solid curve) and ^/s = ITOOGev (dashed curve) 




Q,Gev 

Figure 6: Cronin ratio for McLerran-Venucopalan gluon distribution function for ^4 + ^4 collisions 
for ^/s = 200(761; (solid curve) and i/i = 1700(76^ (dashed curve) 
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Figure 7: Cronin ratio for "dipole" gluon distribution function for A + A collisions for ^/s = 
200Gez; (solid curve) and ^/s = ITOOGef (dashed curve) 





Figure 8: Dependence of ln{qc){foT Raa) on ln{y/s) for different models: Kharzeev-Levin- 
Nardi(solid curve), McLerran-Venugopalan(daslied curve), "dipole" (dot-dashed curve) 
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